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We study Markov exclusion process for a particle system with a local interaction 
in the integer strip. This process models the exchange of velocities and particle- 
hole exchange of the liquid molecules. It is shown that the mean velocity profile 
corresponds to the behaviour which is characteristic for incompressible viscous liquid. 
We prove the existence of phase transition between laminar and turbulent profiles. 
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^ ■ 1 Introduction 

q 

From the mathematical point of view this paper considers a phase transition for a Markov 
multi-dimensional exclusion process. But as this problem originated from concrete physical 
models, we discuss these connections in detail. 

One can study the liquid flow from macroscopic and from microscopic point of view. 
Although in the first case (that is in the continuum mechanics) most studies are based on 
the Navier-Stokes equations, there exist various generalizations introducing randomness in 
these equations (see the review of J. Marsden and other papers in pQ, some history of 
the question one can find in the recent preprint [1]). It is precisely with Navier-Stokes 
equations the immense number of concrete theoretical and practical problems of the liquid 
mechanics, including the Couette flow, which is the simplest (but far from being trivial 
model, see [S]) available for a mathematical study. 

However, there are only very few microscopic models for concrete situations with 
nonzero viscosity, in despite of the abundance of lattice models with infinite number of par- 
ticles (see books [TO], [3], [H])- In theoretical physics the preferable direction is the deriva- 
tion of the Navier-Stokes equations following the sequence: hamiltonian microdynamics 
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— > BBGKY equations — > Boltzmann equation — > equation of hydrodynamics. However, 
there is also direct derivation of the Navier-Stokes equations from stochastic dynamics of 
particles on the lattice (see [XT]). 

Recently, there appeared several papers where some particle system were related to the 
Navier-Stokes equations so that for a certain scaling the dynamics of these particle system 
converges to these equations (see [3], [13], [H]). However, in these papers the particle 
dynamics is not local, it is close to the mean-field dynamics of McKean-Vlasov type. But 
most important is that it is derived not from proper molecular dynamics but from the 
computer simulation of Navier-Stokes equations, that is from the Navier-Stokes equations 
themselves. 

The conclusion is that the derivation of Navier-Stokes equations from physically grounded 
molecular dynamics is still open. The main problem here is that it is far from clear how 
to describe the molecular dynamics. The model of hard balls, obeying the laws of classical 
mechanics, is simple in formulation but very difficult to study. Moreover, one cannot be 
sure that this model is adequate. The question is that the molecules in a liquid are very 
close to each other, and quantum effects seem to be important. If one will work in the 
framework of classical models, then one should take into account the interaction of rotat- 
ing and oscillating degrees of freedom of the neighbor molecules, which presumably have 
a random distributions. 

From the other side, before any attempt to derive the Navier-Stokes equations, it could 
be natural first to try particular cases. In the local model of the Couette flow considered 
below we use two physical processes: exchange of velocities between molecules and the 
exchange between molecules and holes (that is empty places). This corresponds to the in- 
tuitive picture often discussed in the physical literature. Moreover, one can hope that there 
is a kind of universality, that is independence of qualitative behaviour of the concretization 
of the model. 

From probabilistic point of view we study a Markov system of locally interacting par- 
ticles, the dynamics of which is a mixture of two well-known (see [5]) exclusion processes 
- symmetric and completely asymmetric. Introducing boundary conditions does not influ- 
ence the moment closeness of the process [12], and this simplifies essentially the study. 
Precisely using moment closeness one can avoid such difficult techniques as Bethe [TS] or 
matrix [16J ansatz. 

It is worth notice that, due to different scaling, direct study of concrete micromodel 
is not equivalent to the derivation of the Navier-Stokes equations and their subsequent 
application to the same situation. 

For our problem we get the phase transition for the velocity profile, corresponding to 
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the transition in the Couette flow [7j, however there is a difference. The transition pa- 
rameter in our paper is not the classical Reynolds number, which depends on the viscosity, 
cross-section and longitudinal velocity, but its analog, depending on the cross-exchange 
parameter (analog of viscosity), cross-section and random perturbation. The latter could 
be the consequence of fast vortex, resulting from deviations of velocity from longitudinal 
direction or from interaction of velocity with randomly distributed internal degrees of free- 
dom. It appears interesting that the phase transition is sharp on the rougher scale and 
smooth on the finer scale. 

2 Model and results 

The Couette flow is the liquid flow in one (horizontal) direction in the space (which is dis- 
crete and two-dimensional, however the generalization to mult i- dimensional case is straight- 
forward) between two plates. One of the plates is fixed and the other moves with constant 
velocity, that draws the neighbor particles with it. Other particles are also involved due to 
viscosity, that is reflected microscopically in the exchange of velocities between neighbor 
particles. 

We come now to exact definitions. The basic set for us is the discrete strip 

L 8 = {0,l,2,...,S,S + l}xZ. 

To each point (s,x) G L$ we assign random variable vr StX \, which can take one of the 
following three values: 0, or V. The value means that there is no particle at the 
site (s,x) ("a hole"), U( s ,a) — means the presence of one particle with zero velocity, and 
v (s,x) — V means the presence of one particle with velocity V. For notational convenience 
the "hole" will be called the particle with velocity 0. 

Dynamics of this particle system consists of jumps and is defined by the following 
transitions. On any small time interval [t, t + dt] the following events can take place 
independently of one another: 

(a) particles of each neighboring vertical pair (at points (s,x) and (s + l,x)) exchange 
their velocities with probability Xdt + o(dt) (velocity exchange between vertical layers); 

(b) particles in each neighboring horizontal pair of the type V( S)X ) = V, v ( s ,x+i) = ( an d 
only of this type) exchange their velocities with probability \\dt + o (dt) (horizontal flow). 
Of course, one could put V = \i, but in this paper it will be convenient to distinguish V 
and Ai. 

(c) on the zero layer each particle that has velocity V acquires velocity with proba- 
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bility (3 dt + o (dt); on the top layer S+1 each particle that has velocity acquires velocity 
V with probability dt + o (dt) (influence of the boundaries); 

(d) each node (s, x) with nonempty velocity (v( SjX ) = or V) changes its velocity 
according to the rule -H- V with probability edt + o (dt) (random perturbation) . 

In the sequel we assume that A > 0, Ai > 0, (3 > and s > 0. 

Thus we have just defined a continuous time Markov process v(t) = {v( s>x )(t), (s, x) G 
Ls}, t G R+, with state space W = {0, 0, V} Ls . Its generator acts on the functions / that 
depend on velocities at finite number of nodes as follows 



= ) 



(1) 

where two configurations w = v y ^ z and v from W differ only by the velocities at the points 
y and z of Ls' 

Wy = V z , W z = Vy, 

and v y G W is a configuration that can differ from v only in velocity value at node y due 
to the exchange <-¥ V: 

z^y, 

z = y, v y = 0, 

z = y, Vy^0. 

One can check that such definition of C leads to a correct definition of the Markov process 
on W (see Ch. I, §3]). 

Denote by rf the distribution of the stochastic process v(t) at time t. 

Let A = ((si, x{), • • • , (s m , x m )) be some ordered finite subset of Ls and E = (e±, . . . , e m ) 
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be some ordered collection of tj G {0,0, V}. Consider the following probabilities 
q t [A- E) := P {v {suXl) (t) = e 1; . . . , v {Sm>Xm) = e m }. 
On the set L$ we define the action of the group (T a , a G Z) of horizontal translations: 

T a (s,x) = (s,x + a). 

In the present paper we restrict ourself to consideration of initial distributions rf of the 
process v that are invariant with respect to the action of the group T a . 

Since the Markov semigroup generated by ([I]) commutes with the translations T a , for 
any time t the distribution 77* is translation invariant. In particular, for any choice of 
sets A and E 

q t [T a A;E} = q t [A;E] (2) 

for all t > and a G Z. From physical point of view it means that we are interested 
only in homogeneous flows that are invariant with respect to shifts along the x-direction: 
x — > x + a. 

As it is easy to see, the dynamics of the process v is such that the total number of 
particles, namely, the number of nodes with velocities and V, is a conserved quantity. 
In particular, the mean (expected) number of nonempty nodes in a vertical section x is 
conserved in time and does not depend on x: 

s+i 

J2(Qt[{k,x);0} + q t [{k,x);V}) = M Vx,t. (3) 

Evidently, the constant M = M(rf) can be calculated in terms of the initial distribution rf. 
However, the mean number of nodes with velocity V is not conserved. 

From now on we are interested only in distributions in a fixed vertical layer, therefore, 
via the homogeneity, the dependence on x will be omitted in the following notation 

q t [(s,x);e]=pt(t). (4) 

Theorem 1 Assume that the initial distribution rf is translation invariant. Then: 
1) {existence of a stationary regime) the following limits exist for any s and e 



\i e s = lim p e s (t), 

t — ^00 
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we call [L e s stationary probabilities; 

2) (uniformity of a particle density) the stationary probabilities of "holes" are the same 
for any layer s = 0, 1, . . . , S + 1: 

/if = const. 

First we consider the case when the random perturbation is absent, i.e., e = 0. 

Theorem 2 (linear velocity profile) Assume that the initial distribution rj° of the process 
is translation invariant and e = 0. Then the mean velocities of particles form a linear 
profile, namely, 

^ = g + 1 + 2A/3- 1 (fc + A/rl) ' k = 0,l,...,S + l, 

where p s = M(rf)/(S + 2). 

The velocity profile of this theorem corresponds to the laminar flow in contrast to what 
is observed in the turbulent regime (see [7]). In the next theorem we consider nonzero 
random perturbations. 

Theorem 3 Assume that S — > oo, the parameter e = e$ depends on S and the sequence 
of initial distributions {i] 0,s } is such that the density of particles is fixed, i.e., 

M^^p 6 [0,l], S^oo. (5) 

Then: 

1) (rough scale picture of the transition) in scaling s = uS the functions //° and p7 s 
demonstrate a phase transition from laminar profile to turbulent one; namely, ifesS 2 — > 0, 
then the limiting profile of the function fj^ s , u G (0, 1), has a linear (laminar) form: 

vis -> P u ( s -> °°). 
but, if EsS 2 — > oo, then the limiting profile is constant (turbulent): 

Pus |> ( s -> oo); 



2) (finer scale picture of the transition) the case EsS 2 = const, S — > oo, is intermediate 
between the situations described in point 1). Precisely, if we introduce a positive parameter 
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K putting esS 2 



^XK 2 , then — y gx(u), where 




1 + 



smh(K(u - 1/2)) 
sinh(AT/2) 



) 



Since \im K ^ g K (u) = pu and hin^^oo gx(u) = p/2, we can say that by changing the 
parameter K > we observe a smooth transition between laminar and turbulent profiles. 

Remark 1. The parameter K can be considered as analog of the Reynolds number 
Re . This can be easily explained by comparing these two numbers, that is by comparing 
the definition of K 



with the definition of the Reynolds number 



where L is the cross-section of the flow, V is the mean velocity, v is the viscosity. 

The mentioned phase transition is based on the competition between processes which 
reinforce and the process which diminish the influence of the boundary Thus, the fac- 
tors L and S show that the greater the distance from the boundaries, the more random 
perturbation diminishes the influence of the boundaries. The factors V and e show that 
if the velocity increases the particles are perturbed more often, and the influence of the 
boundaries is less visible. Finally, when A increases, then the influence of the boundaries 
increases. One should note that, by definition of the model, the parameter V can be put 
proportional to Ai. 



3.1 Closed equations for one-particle functions 

Our next task is to obtain a system of differential equations for the functions pl{t), < 
k < S + 1, e G {0, 0, V}, defined in (0}. First note that for any k = 0, 1, . . . , S + 1 



Re 



LV 



v 



3 Proofs 



P k (t)+p° k (t)+p v k (t) = l w>o 



(6) 



hence, for example, the functions p° k {t) can be easily expressed in terms of the functions 
pf(t)and p\{t). 
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Lemma 1 Let the initial distribution rf be translation invariant. Then the functions p e k {t) 
satisfy to a closed system of linear first-order differential equations having the following 
structure. 

1) The functions p k (t) can be obtained from a subsystem 

X(p k+1 (t)+ P ti(t)-2p k (t)), l<k<S, 

X(p (t)-p (t)), (7) 
Hp s(t)-P s+I (t)). 

2) The equations for p\ (t) have the following form 

Kpl + i(t)+pU(t) - 2p v k (t))+e(l-p k (t) - 2pX(t)), 
l<k<S } 

A(pT(f) - P V a (*)) + e(l - p (t) - 2p%(t)) - 0pg(t), (8) 

Kp V s(t) -P v s + i(t)) +e(l -p% +l {t) - 2p v s+1 (t)) 
+ P(l-P s+I (t)- P v +I (t)). 

Note that equations (7) and (8) do not depend on Ai. 

Remark 2. It can easily be checked that the system ([7]) has the same form as the corre- 
sponding system for the simple symmetric exclusion process on the finite set 0, 1, . . . , S + 1 
(with empty boundary conditions). 

The remaining part of this subsection is devoted to the proof of Lemma [TJ First we 
shall obtain equations for the marginal distributions q t [(k,x);e\ without assumption of 
translation invariance of the initial distribution of velocities. 



d 
~dt 
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For inner layers 1 < k < S the marginal distributions satisfy the following equations. 

j t q t [(k, x); 0} = X(q t [(k + 1, x); 0] + q t [(k - 1, x); 0] - 2q t [(k, x); 0]) 

+ Ax (q t [((k, x), (k, x + 1)); (V, 0)] - q t [((k, x-1), (k, x)); (V, 0)]) , (9) 

n 

dt 



^ q t [(k, x); 0] = X(q t [(k + 1, x); 0] + g t [(fc - 1, x); 0] - 2q t [(k, x); 0]) 



+ e(<fc[(M);^]-ft[(M);0]), 
^ q t [(k, x); V] = X(q t [(k + 1, x); V] + g t [(fc - 1, x); V] - 2q t [(k, x); V]) 
+ e{q t [(k,x);0]-q t [(k,x);V}) + X, ( - q t [((k, x), (k, x + 1)); (V, 0)] 

+ &[((M-l),(M));G / ,0)]). (10) 

On the bottom layer k — the equations have the following form 

^&[(O,x);0] = X(q t [{l,x)-0\ -q t [(O,x);0]) 

+ X 1 (q t [((0,x), (0,x + 1)); (V, 0)] - %[((0,x - 1), (0,x)); (V, 0)]), (11) 

j t q t [(0, x); 0] = A(g t [(l, x); 0] - g t [(0, x); 0]) + ^ t [(0, x); V] 
+ e(q t [(0,xy,V]-q t [(0,x);0}), 

j t q t [(0, x)-V} = A(g t [(l, x); V] - g t [(0, x); V]) - #ft[(0, x); V] 

+ £(- gt [(0,x);F]+ ?t [(0,x);0])+ A 1 (-gi[((O,x),(O,x + l));(\/,0)] 

+ gt [((O,x-l),(O,x));(y,0)]). (12) 
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On the top layer k = S + 1 we have the equations 



d 
~dt 



q t [{S + l,x);0] 



X(q t [{S,xy,0}-q t [(S + l,x);0}) 

+ \ 1 (q t [((S + l,x),(S + l,x + l));(V,0)\ 



d 
dt 



d 
dt 



q t [(S + l,x);V] 



q t [(S + l,x);0] 



-q t [((S + l,x-l),(S + l,x));(V,0)]), 
X(q t [(S, x); 0] - q t [(S + 1, x); 0]) + (3q t [(S + l,x); V] 
+ e(q t [(S + l,xy,V]-q t [(S + l,xy,0}), 
X(q t [(S, x); V] - q t [(S + l,x); V]) + (3q t [(S + 1, x); 0] 
+ e(q t [(S + l,x);0]-q t [(S + l,x);V]) 
+ A x ( - q t [((S + l,x),(S + l,x + 1)); (V, 0)] 



(13) 



+ q t [((S + l,x-l),(S + l,x));(V,0)]). 



(14) 



We see that the equations for one-dimensional marginal distributions are not closed, 
since they contain two-dimensional distributions in the r.h.s. 

Derivation of these equations is straightforward, but intermediate calculations are 
rather cumbersome. This derivation will be more transparent if we note that dynam- 
ics of the model consists of three processes. The first one is a symmetric exclusion process 
on vertical layers. The second component is a totally asymmetric exclusion process in the 
horizontal direction. And, finally, the Glauber spin-flip process at each point. This latter 
process enters the equations in a very simple manner, but it is very useful to consider the 
first two processes separately. We shall see that the equations for the first process are 
closed even with the boundary conditions, while the equations for the second process are 
not closed. 

Model 1: symmetric exclusion process. First it is convenient to consider a simpler 
auxiliary model £(t) = (£ s (£)), which describes a single vertical layer x = 0, consisting of 
points s — 0, 1, 2, . . . , S, S + 1, under additional assumption Ai = 0, e = and (3 = 0. For 
notational convenience we assume that there are no "holes" , i.e., £, s (t) G {0, 1}, and put 
V = 1. Then the process t > 0, is a well-known simple exclusion process. 

For l<i<i + k — 1 < S and e,j = 0, 1 denote 



One-particle functions pt(s; 1) satisfy closed equations. Indeed, using an explicit form of 
the transition functions of the Markov process £ on a small time interval [t, t + dt) and the 



p t (i; e 1 . . . e k ) = P = e u . . . , &+k_i(*) = e k }. 
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complete probability formula, for any 1 < % < S we have, up to terms of order o (dt), 

p t+dt (i;l) = p t (i-l;lll) + [p t (i-l;011)+p t (i-l;110)](l-AtZt) 
+ p t (i - 1; 010)(1 - 2A dt) + \p t (i - 1; 100) 

+ p t {i - 1; 001)] A dt + p t {i - 1; 101) • 2A dt, (15) 

hence 

Pt+dt(i; 1) - Pt(i] 1) = Adt [-p t (i - 1; 01) - p t (i; 10) + p t (i - 1; 10) + p t (i; 01)] 

= \dt\p t {i-l;l)-2p t (i;l)+p t (i + l;l)}, (16) 

and, finally, 

j t p t (i; 1) = A[p t (i - 1; 1) - 2p t (i; l)+p t (i + 1; 1)]. (17) 

Model 2: asymmetric exclusion process. It is convenient also to consider an asymmetric 
exclusion process ((t) = (( x (t), x e Z) with three values at the node: Cx £ { ,O, V}. We 
assume that the following transitions are only possible: each pair of closest neighbors of 
the type V0 become 0V with intensity Ai. In other words, on the time interval [t, t + dt] 
each particle, having velocity V, independently of other particles jumps one unit to the 
right (x — > x+ 1) with probability Ai<it + o (dt) provided that the node x + 1 is empty. We 
need to do calculations similar to (TToT) and (TTBT) . Now x G Z. As before, let us consider 
all possible states in the neighboring nodes x — 1, x and x + 1 at time t and apply the 
complete probability formula: 

p t+dt (x; 0) = J2 P M* + ^) = I (Cx-i(t), CxW, C*+i(*)) = (ei, e 2 , e 3 )) 

eie{0,O,y},j=l,2,3 

xP{(C x ._i(t),CxW,Cx+iW) = (ei,e 2 ,e 3 )}. 

Conditional probabilities under the sum sign can be easily found. The results are presented 
in the following table, where "*" means any of symbols 0, or V: 



x-1 


X 


x + 1 


P(Cx(t + dt) = £ 


3|(C«-i(*).C*(t),C«+i(*)) = (ei,e a ,e 3 )) 


ei 


e2 














* 




1 - o (<it) 












1 - o (dt) 


V 









1 - Ai dt + o (dt) 


* 







o(dt) 


* 


V 







X 1 dt + o (dt) 


* 


V 







o(dt) 


* 


V 


V 




o(dt) 
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Up to terms of order o (dt) we get 



Pt(x - 1; 00) +p t {x- 1; 00) + pt(aj - 1; V0) 
- Ai dt ■ p t (x - 1; 1/0) + Ai ■ pt(a;; V0) 
Pt(ar; 0) + A x dt • [-p t (x - 1; V0) + pt(ar; V0)] 




(18) 



d 
dt 



pt(x; V) = Ai[pt(z - 1; V0) - p t (x; F0)]. 



Similarly one can consider the components of the dynamics (pQ), corresponding to the 
random perturbations and to the behavior on the boundaries. 

It appears that the equations for q t [(k,x);e] become closed if we consider only trans- 
lation invariant distributions. Namely, to finish the proof of Lemma HI from now on we 
assume that the process starts at time t = from a translation invariant distribution. 
Recall that the dynamics conserves the translation invariance, i.e. at any time t > 
the process has a translation invariant distribution. In particular, the property (j2J) holds. 
Therefore, in the right hand sides of the equations (j91)- (IT3|) all summands having the factor 
Ai disappear. Collecting together the equations with q t ((k,x); 0), < k < S + 1, and 
using notation (jlj) , we come to the statement 1) of Lemma HJ The statement 2) easily 
follows from the form of the equations (151)- (1141) if we take into account ([H]). 

This completes the proof of the lemma. 

3.2 Convergence as t — > oo 

We prove here case 1) of Theorem Uj namely, we show the existence of the limits lim^oo p%(t) 
for the functions p%(t), satisfying the equations (JTJ)— (|HD - We use probabilistic arguments. 
On the finite state space M. = {0, 0, y r }{ > 1 >--»S'+i} us define an auxiliary continuous 
time Markov process 



= (#*(*), k = 0,l,...,S + l) 



with the following transitions. On a small time interval [t, t + dt] the following events can 
occur independently of each other: 
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- any nodes k and k + 1 exchange their states with probability A dt + o (dt); 

- any node of the zero layer k = 0, being in the state V, changes its state to with 
probability (3 dt + o (dt); any node of the top layer k = S + 1, being in the state 0, with 
probability (3 dt + o (dt) changes its state to V; 

- any node k in a state different from changes its state according to the rule •<-)■ V 
with probability edt + o (dt) (random perturbation) . 

Note that the Markov process $(t) is reducible, because its dynamics does not change 
the total number of nodes having the empty state (0). At the same time the process i)(t) 
considered on any set, invariant with respect to the dynamics 

M n = {&=(# ,...,>&s+i)eM: # j = 0} = n}, n = 0, 1, . . . ,3 + 1, 

is irreducible and ergodic. Hence, for any n the following limits 

lim P (d(t) = x | t?(0) = x ) = g n (x), xo, x G M n , 

exist and do not depend on the concrete choice of initial state x . Thus for any initial 
distribution /z#(o) there exists the limit lim^oo P {$(t) = x}, depending on (i$(o). 
If we permit, within the current subsection, the following notation 

p e k (t) = P {& k (t) = e}, ee {0, 0, V}, k = 0, 1, . . . , S + 1, 

then it is easy to see that the functions p%(t) satisfy the system of equations ©-([8]). Proof 
of this fact is even simpler than the reasonings of subsection 13. 11 Since 

Pt(t)= £ p{*(f) = x}, 

the statement 1) of Theorem [1] easily follows. 

3.3 Exact formulas for stationary solutions 

Now we return to the process v(t). Denote by /if, u° s and the stationary probabilities of 
the following three events: on the vertical layer s there is no particle, there is a particle with 
velocity and there is a particle with velocity V correspondingly, that is lim p e s (t) — n e s . As 
it was explained in the previous subsection these limits depend on the initial distribution rf. 
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In particular, 

s+i 

= + 2) -Mfo°), (19) 

s=0 

where M(rj°) is defined in ([3]). 

Let us remark that the limiting probabilities pf s satisfy the stationary versions of the 
equations (j7|), ([S]). It follows from (JB} that 

tf + IM°, + nY = l Vs. 

Lemma 2 T/ie distribution of "holes" is uniform on the set of layers: /if = /if /or all 

s, Sl e {0,1,..., 5 + 1}. 

Proof The system of stationary equations for /if , s = 0, 1, . . . , S + 1, has the form: 

= A&i+A&i-Spf, l<s<5, 
= /if -/if, 

It is evident that the set of all solutions of this system is the one- dimensional subspace 
/if = • • • = /if, r From the probabilistic nature of our model we have 

//f = ... = pf +1 = l-p s e[0,i]. 

This proves the lemma and the statement 2) of Theorem [TJ 

The number p$ G [0, 1] will be called particle density. It is readily seen from f fl9|) that 
the density depends on the initial distribution: p s = M(rf)/{S + 2). 

Since /i° + /i^f = ps, we need only to find probabilities (/ijf, fc = 0, 1, . . . , S + 1). On 
the inner layers they satisfy the following equations: 

= \(i4 +1 + t4_ 1 -2^)+e(p s -2f4), l<k<S, (20) 

with the following boundary conditions 

= -^)+e( Ps -2^)-^, (21) 

= X(p v s -p v s+l )+6(ps-2p^ +1 ) + P(p s -p v s+l ). (22) 



Case 1: e — 0. Equations (1201 take the following form /i^ +1 + p^-i ~ = 0. 
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Characteristic equation corresponding to this difference equation 

z 2 + 1 - 2z = 

has a root z — 1 of order 2. This is the resonant case, therefore a general solution is 
= Do+Dik. The coefficients Do and .Di can be found from the boundary conditions ( I2ip 
and (F22]): 

= XDi - 0D O , = -XD 1 + f3( Ps -D -D 1 (S+ 1)). 
Solving the latter system, we get explicit solution 

^ = ^^xpi^ + ^' 1 ^ k = 0,l,...,S + l, 

which is a linear function. Now the statement of Theorem |2] easily follows. 

Case 2: e > 0. It can be shown in a standard way that general solution of the inhomo- 
geneous difference equation ( f20l) has the form 

fik = y + C x z\ + C 2 z k 2 , (23) 
where z\ and zi are the roots of the characteristic equation 

A^ 2 + A-2(A + e)^ = 0, (24) 
C\ and C*2 are unknown coefficients. Immediately, by solving equation (|24|) . we get 



Zi, 2 = l + |±A/fl + | N ) 2 -l. (25) 



A y v A 



Notice that if e > then the equation (1241) has two different real roots: zi > 1 and z% < 1, 
moreover, = 1- Substituting (123]) to the boundary conditions (pT!) - (!22|) . after some 
algebra we obtain the next lemma. 



Lemma 3 The solution (/i^, s = 0,l,...,S' + l) has the following explicit f 



orm: 



v p s (&(^f +1 ~* - b{z 2 )z s 2 +1 - k ) - (b( Zl )z k - b(z 2 )z k 2 ) 



2 ^ 2(6 2 (z 1 )zf +1 -6 2 (z 2 )z 2 5+1 ) 
fc = 0, 1, . . . , 3 + 1, where b(z) := 1 + (2e + A(l - z" 1 ))//?. 
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3.4 Asymptotics and phase transition 

Here we consider A and (3 as fixed parameters and assume that e —> 0, S — > oo. We see 
that for small e 

Taking into account the assumptions of Theorem [3j for simplicity we put ps = p- 
We are interested in the behaviour of , < k < S + 1) on the space scale k = [uS], 

u E [0,1]. To get a meaningful asymptotics in the limit S — > oo, we assume that the 

parameter e is a function of S: 

es = ^\S)S->, 

where ijj(S) > is such that ip(S)S~ 1 — > as S — > oo. Hence, yj2es/\ = ip(S)S~ 1 and 

z 1)2 = 1±^(S , )S , - 1 + 0(^ 2 (S , )S'- 2 ), 

b(z 1>2 ) = i±^{s)s- 1 + r 1 o{^ 2 {s)s- 2 ). 

The calculations presented below show that the asymptotics of pY u s\ strongly depends on 
the properties of the function ip. 

In case of small random perturbations ip{S) — > 

»W ~ 2- p 2~2W) + /? ° {S ] 

p 2u — 1 
~ 2 +P — 2 — =MP ' 

and the profile has a laminar character. 

In case of strong random perturbations ip(S) — > oo 

(l-w)S _ u s 

?U " 2 2z s " 2 + r l O{S-" ■ ^S)e~^) WE (0, 1). (26) 

and the limiting profile is very far from being laminar. Indeed, it can be checked that 
the limiting profile here is the same as in the case when there is no boundary conditions 
(formally, when (3 = 0), that is called free boundary conditions. We can conclude that if 
the random perturbations are strong then the boundaries do not influence the inner layers. 
This is typical for the turbulent phase. 
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The case 4>(S) = const separates the above two phases. For convenience let us put 
tp(S) = K,K>0. Then 

v p (e K ^- u ) - e -^( 1 -")) - ( e ^ - e~ Ku ) 



= P/ 1 + sinh(i,(,-l/2)) x 
2v + sinh(A72) /' 

Let us remark in conclusion that the assumption -^(S^S 1-1 — > is not essential for the 
strong perturbation asymptotics fl26|) . In particular, this results holds in the important 
case es = £■ Indeed, in this case z 2 does not depend on S, and, using the explicit formula 
of Lemma [31 we get for any u G (0, 1) 

v P const ^ n (Q _± ^ 

^[ui>\ q r 7 / x (i_ m m(l— u,u))S 
2 



b(z 2 )z 



We stress that all limiting expressions do not depend on the parameter /3. 
Theorem [3] is proved. 



References 

[1] Strange attractors. Mir Moscow, 1981. (Mathematics. News of foreign science, V. 22.) 

[2] Chorin A. J., Marsden J. E. A Mathematical Introduction to Fluid Mechanics. 
New York: Springer- Verlag, 1993, 169 p. 

[3] Marchioro C, Pulvirenti M. Mathematical Theory of Incompressible Non Vis- 
cous Fluids. New York: Springer- Verlag. 1994, 283 p. 

[4] Shiryaev A. On the classical, statistical, and stochastic approaches to the hydrody- 
namic turbulence. Research Report N.2. Aarhus: Thiele Centre for Applied Mathe- 
matics in Natural Science, 2007. 

[5] Kolmogorov A.N. Mathematical models of turbulent motion of an incompressible 
viscous fluid. RUSS MATH SURV, 2004, v. 59, N.l, p. 3-9. 

[6] Meshalkin L. D., and Sinai Ia. G. Investigation of the stability of a stationary 
solution of a system of equations for the plane movement of an incompressible viscous 
liquid. Journal of Applied Mathematics and Mechanics, 1961, v. 25, N.6, p. 1700-1705 



18 Stochastic Micromodel Of The Couette Flow 

[7] Landau L.D., Akhiezer A. I., Lifshitz E.M. General Physics: Mechanic and 
Molecular Physics. Pergamon Press, 1970, 372 p. 

[8] Liggett T. Interacting Particle Systems. Springer, 2005, 496 p. 

[9] Kipnis C, Landim C. Scaling Limits of Interacting Particle Systems. Berlin: 
Springer- Verlag, 1999, 442 p. 

[10] Spohn H. Large Scale Dynamics of Interacting Particles. Berlin: Springer- Verlag. 
1991, 342 p. 

[11] Esposito R., Marra R., Yau H. T. Navier-Stokes equations for stochastic particle 
systems on the lattice. Comm. Math. Phys., 1996, v. 182, N.2, p. 395-456. 

[12] Ignatyuk I.A., Malyshev V. A., Molchanov S.A. Moment-closed processes 
with local interaction. Selecta Math. Soviet., 1989, v. 8, N.4, p. 351-384. 

[13] Meleard S. Monte- Carlo approximations for 2d Navier-Stokes equations with mea- 
sure initial data. Probab. Theory Related Fields, 2001, v. 121, N.3, p. 367-388. 

[14] Fontbona J. A probabilistic interpretation and stochastic particle approximations 
of the 3-dimensional Navier-Stokes equations. Probab. Theory Related Fields, 2006, 
v. 136, N.l, p. 102-156. 

[15] Golinelli O., Mallick K. Derivation of a matrix product representation for the 
asymmetric exclusion process from algebraic Bethe ansatz. J. Phys. A., 2006, v. 39, 
N.34, p. 10647-10658. 

[16] Derrida B., Lebowitz J.L., Speer E. R. Entropy of open lattice systems. 
J. Statist. Phys., 2006, v. 126, N.4/5, p. 1083-1108. 



